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Abstract 

A general method is known to exist for studying Abelian and non-Abelian gauge theories, as well 
as Euclidean quantum gravity, at one-loop level on manifolds with boundary. In the latter case, 
boundary conditions on metric perturbations h can be chosen to be completely invariant under 
infinitesimal diffeomorphisms, to preserve the invariance group of the theory and BRST symmetry. 
In the de Bonder gauge, however, the resulting boundary-value problem for the Laplace type 
operator acting on h is known to be self-adjoint but not strongly elliptic. The latter is a technical 
condition ensuring that a unique smooth solution of the boundary-value problem exists, which 
implies, in turn, that the global heat-kernel asymptotics yielding one-loop divergences and one- 
loop effective action actually exists. The present paper shows that, on the Euclidean four-ball, 
only the scalar part of perturbative modes for quantum gravity are affected by the lack of strong 
ellipticity. Further evidence for lack of strong ellipticity, from an analytic point of view, is therefore 
obtained. Interestingly, three sectors of the scalar-perturbation problem remain elliptic, while lack 
of strong ellipticity is "confined" to the remaining fourth sector. The integral representation of the 
resulting ^-function asymptotics is also obtained; this remains regular at the origin by virtue of a 
spectral identity here obtained for the first time. 

PACS: 03.70.+k, 04.60.Ds 
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I. INTRODUCTION 



The study of gauge theories and quantum gravity on manifolds with boundary is mo- 
tivated by the problems of quantum cosmology il| and quantum field theory under the 
influence of external conditions 0], and several efforts have been produced in the literature 
;o study boundary conditions and one-loop semiclassical properties within this framework 
^. In our paper we are interested in boundary conditions for metric perturbations that are 
completely invariant under infinitesimal diffeomorphisms, since they are part of the general 
scheme according to which the boundary conditions are preserved under the action of the 
symmetry group of the theory . In field-theoretical language, this means setting to 

zero at the boundary that part ttA of the gauge field A that lives on the boundary B (tt 
being a projection operator): 



as well as the gauge-fixing functional, 



and the whole ghost field 



ttA 



0, 



0. 



;i.2) 



:i.3) 



For Euclidean quantum gravity, Eq. (1.1) reads as 

= 0, (1.4) 

where hij are perturbations of the induced three-metric. To arrive at the gravitational 
counterpart of Eqs. (1.2) and (1.3), note first that, under infinitesimal diffeomorphisms, 
metric perturbations h^^i, transform according to 

where V is the Levi-Civita connection on the background four-geometry with metric g, 
and Lp^dx^ is the ghost one-form (strictly, our presentation is simplified: there are two 
independent ghost fields obeying Fermi statistics, and we will eventually multiply by —2 the 
effect of (fij to take this into account). In geometric language, the infinitesimal variation 



5h 



h 



hfj^u is given by the Lie derivative along (p of the four-metric g. For manifolds 



with boundary, Eq. (1.5) implies that 



m 



hij = hij + ip(^i\j) + Kijipo, 



[1.6) 



where the stroke denotes three-dimensional covariant differentiation tangentially with re- 
spect to the intrinsic Levi-Civita connection of the boundary, while Kij is the extrinsic- 
curvature tensor of the boundary. Of course, (/?o and are the normal and tangential 
components of the ghost, respectively. By virtue of Eq. (1.6), the boundary conditions (1.4) 
are "gauge invariant", i.e. 

%1 =0, (1.7) 
if and only if the whole ghost field obeys homogeneous Dirichlet conditions, so that 



OB 



0, 



:i.8) 



Nb = 0. (1.9) 

The conditions (1.8) and (1.9) are necessary and sufficient since lpq and (fi are independent, 
and three-dimensional covariant differentiation commutes with the operation of restriction 
to the boundary. We are indeed assuming that the boundary B is smooth and not totally 
geodesic, i.e. Kij ^ 0. However, for totally geodesic manifolds, having Kij = 0, the condition 
(1.8) is no longer necessary. 

On imposing boundary conditions on the remaining set of metric perturbations, the key 
point is to make sure that the invariance of such boundary conditions under the infinitesimal 
transformations (1.5) is again guaranteed by (1.8) and (1.9), since otherwise one would obtain 
incompatible sets of boundary conditions on the ghost field. Indeed, on using the DeWitt- 
Faddeev-Popov formalism for the (out | in) amplitudes of quantum_gravity, it is necessary to 
use a gauge-averaging term in the Euclidean action, of the form [2| 

1 r <^,.<i>^ 



where is any functional which leads to self-adjoint (elliptic) operators on metric and 
ghost perturbations. One then finds that 

5<^,ih) ^ <^,ih) - <!>,(h) = ^; v^,, (1.11) 

where J-'^^ is an elliptic operator that acts linearly on the ghost field. Thus, if one imposes 
the boundary conditions 
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= 0, (1.12 

J B 



and if one assumes that the ghost field can be expanded in a complete orthonormal set of 
eigenfunctions u^^^ of J^^ which vanish at the boundary, i.e. 

^>i'^ = An(^\ (1.13) 



(A) 



;i.i4) 



u 



(A) 



0, 



;i.i5) 



the boundary conditions (1.12) are automatically gauge-invariant under the Dirichlet con- 
ditions (1.8) and (1.9) on the ghost. 

Having obtained the general recipe expressed by Eqs. (1.4) and (1.12), we can recall what 
they imply on the Euclidean four-ball. This background is relevant for one-loop quantum 
cosmology in the limit of small three-geometry on the one hand Q|, and for spectral geom- 



etry and spectral asymptotics on the other hand 
the de Bonder gauge-fixing functional 



1 Ij . As shown in Ref. 



0, 



if one chooses 



:i.i6) 



which has the virtue of leading to an operator of Laplace type on h^j^y in the one- loop 
functional integral, Eq. (1.12) yields the mixed boundary conditions 



^ + 7^oo-^(<7^/ii,) + ;^/i(j. 



0, 



:i.i7) 



J B 



dhpi _^ 3^ _ 1 dhpo 
dr r °* 2 dx'^ 



0. 



:i.i8) 



In Refs. j3|,|7|], the boundary conditions (1.4), (1.17) and (1.18) were used to evaluate the 
full one-loop divergence of quantized general relativity on the Euclidean four-ball, including 
all h^jj and all ghost modes. However, the meaning of such a calculation became unclear 
after the discovery in Ref. |6j that the boundary-value problem for the Laplacian P acting 
on metric perturbations is not strongly elliptic by virtue of tangential derivatives in the 

n 

boundary conditions (1.17) and (1.18). Strong ellipticity [11] is a technical requirement 
ensuring that a unique smooth solution of the boundary-value problem exists which vanishes 
at infinite geodesic distance from the boundary (see Appendix A). If it is fulfilled, this ensures 
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that the LP' trace of the heat semigroup e~*^ exists, with the associated global heat-kernel 
asymptotics that yields one-loop divergence and one-loop effective action. However, when 
strong ellipticity does not hold, the trace of e~*^ acquires a sing 
it is unclear how to attach a meaning to ^-function calculations. 

All of this has motivated our analysis, which therefore starts in Sec. II with the mode- 
by-mode form of the boundary conditions (1.4), (1.8), (1.9), (1.17) and (1.18) with the 
resulting eigenvalue conditions. Section III studies the matrix for coupled scalar modes, 
while Sec. IV obtains the first pair of resulting scalar-mode C-functions and Sec. V studies 
the remaining elliptic and non-elliptic parts of spectral asymptotics. Results and open 
problems are described in Sec. VI, while technical details are given in the Appendices. 



II. EIGENVALUE CONDITIONS ON THE FOUR-BALL 



On the Euclidean four-ball, which can be viewed as the portion of fiat Euclidean four- 
.pace bounded by a th.ee-sphere of .^up^ mat* perturbations can be expanded in 
terms of hyperspherical harmonics as 



aams q, 

oo 

/ioo(x,r) = ^a„(r)Q(")(x), 



n=l 



n=2 



bn{r) 



[n^ — 1) 



(2.1) 



(2.2) 



hij{x,T) = ^(i„(r) 



n=3 
oo 



X] 



oo , s 



n=l 



n=3 



+ E \fni'r)(sl;^{x) + 5];)(x) + K{t)G[ 



XI 



(2.3) 



where r G [0, q] and g(")(x), 5f ^(x) and Cj"^ (x) are scalar, transverse vector and transverse- 
traceless tensor hyperspherical harmonics, respectively, on a unit three-sphere with metric 
Cij. By insertion of the expansions (2.1)-(2.3) into the eigenvalue equation for the Laplacian 
acting on /i^^,, and by setting \/E — > iM, which corresponds to a rotation of contour in 
the C-function analysis 141, one finds the modes as linear combinations of modified Bessel 
functions of first kind according to 



a„(r) = - 7i/„(Mr) +73/„_2(Mr) +74J„+2(Mr) 

T I 



(2.4) 



bn{r) = 72/„(Mr) + {n + l)73/„_2(Mr) - (n - l)74/„+2(Mr), 



(2.5) 



c„(r) = (Mr) + e^In-iiMr) 



(2.6) 



72/n(Mr) + i!^^3J„_2(Mr) + i!^74/„+2(Mr) 



(n-2) 



in + 2) 



(2.7) 



e„(r) = r (871 - 272)4 (Mr) - 734_2 (Mr) - 744+2 (Mr) 



(2.8) 



/n(r) = r 



£1 



£2 



(r2 + 2) 



4+1 (Mr) + — ^4_i(Mr 
{n-2) 



(2.9) 



kn{T)=aiTln{MT). (2.10) 

Modified Bessel functions of second kind are not included to ensure regularity at the origin 
r = 0. Moreover, normal and tangential components of the ghost field admit the following 
expansion on the four-ball: 



V.o(x,r) = ^/„(r)Q(")(x), 



n=l 



(^i(x,r) = ^ 



n=2 



^n[r)— ^^Vn[j)S\ \X) 



where the ghost modes /„(''"), fnJyT^ and p„(r) are found to read as 



ln{T) = - Ki4+i(Mr) + K24-i(Mr) | , 



(2.11) 



(2.12) 



(2.13) 



m„(r) = -{n - l)«:i4+i(Mr) + (n + l)«:24-i(Mr), 



(2.14) 



Pn{r)=ein{MT). (2.15) 

At this stage, the boundary conditions (1.4), (1.17), (1.18), (1.8) and (1.9) can be re- 
expressed in terms of metric and ghost modes as 



dan Q I dcn 2 3 

-} h -an 5--} = on 5 , 

ar T ar r^ 
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(2.16) 



dbn 3 , 
-r- + -On 
ar T 



in' - 1: 



■a„ = on 5^ 



(2.17) 



,3 3 

—, 1 — c„ = U on 6 , 

ar r 



(2.18) 



dn = en = fn = K = L = mn = Pn = on S"^. (2.19) 

Furthermore, the formulae (2.4)-(2.10) and (2.13)-(2.15) can be used to obtain homogeneous 
hnear systems that yield, implicitly, the eigenvalues of our problem. The conditions for find- 
ing non-trivial solutions of such linear systems are given by the vanishing of the associated 
determinants; these yield the eigenvalue conditions S{E) = 0, i.e. the equations obeyed by 
the eigenvalues by virtue of the boundary conditions. For the purpose of a rigorous analysis, 
we need the full expression of such eigenvalue conditions for each set of coupled modes. 
Upon setting ^/E — > iM, we denote by D(Mq) the counterpart of S{E), bearing in mind 
that, strictly, only d{E) yields implicitly the eigenvalues, while D{Mq) is more convenient 



for ^-function calculations 



To begin, the decoupled vector mode C2(r) = I^{Mt) obeys the Robin boundary condition 
(2.18), which yields 

D{Mq) = h{Mq) + h{Mq) + -^h{Mq), (2.20) 

Mq 

with degeneracy 6. Coupled vector modes c„(t) and fnij) obey the boundary conditions 
(2.18) and (2.19), and hence the corresponding D{Mq) reads as 

Dn{Mq) = 4_i(Mg)(j„(Mg) + 4+2(Mg) + ^J„+i(Mg)) 

+ l^rf (4-2(Mg) + 4(Mg) + A4_^(Mg)) , (2.21) 

with degeneracy 2{n'^ — 1), for all n > 3. 
The scalar modes 

ai(r) = i(7i/i(Mr) +74/3(Mr)), (2.22) 

ei(T) = r(37i/i(Mr) -74/3(Mr)), (2.23) 
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obey the boundary conditions 



dai 6 1 dei „ ,^ 
+ ^-7^ = at T = g, 2.24 



ei(g)-0, 



(2.25) 



which imply 



D{Mq) = 2Qh{Mq)h{Mq) - Mq{Io{Mq) + h{Mq))h{Mq) 
+ 3Mqh{Mq){l2{Mq) + h{Mq)), 



with degeneracy 1. 
The scalar modes 



a2 r 



7i/2(MT)+74/4(Mr) 



(2.26) 



(2.27) 



62(r) =72/2 (Mr) -74/4(Mr), 



(2.28) 



e2(T) = r (371 - 272)/2(Mr) - -fMMr) 



(2.29) 



obey the boundary conditions (2.16), (2.17) and (2.19) with n = 2, and hence yield the 
determinant 



D{Mq) = det 



hiMq) - Mql!^{Mq) Mql!^{Mq) AhiMq) + MqIi{Mq) 

ShiMq) -{2MqI^{Mq) + Qh^Mq)) 2Mqr^{Mq) + 9/4(Mg) 

ShiMq) -2h{Mq) -h{Mq) 

(2.30) 
with degeneracy 4. 

For all n > 3, coupled scalar modes a„,6„,(in,e„ obey the boundary conditions (2.16), 
(2.17), (2.19). The resulting determinant reads as 



Dn{Mq)^dQi pij{Mq), 
with degeneracy n^, where pij is a 4 x 4 matrix with entries 



(2.31) 



pn = In{Mq) - Mql'^iMq), pu = MqI'^{Mq), 



(2.32) 
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Pi3 = (2 - n)In-2{Mq) + MqI'^_^{Mq), pu = (2 + n)4+2(Mg) + MqI^_^^{Mq), (2.33) 



P21 = -(n' - l)/n(Mg), P22 = 2MqI'^{Mq) + 64(Mg), (2.34) 



P23 = 2(n + l)Mgj;_2(Mg) - (n^ - 6n - 7)4_2(Mg), (2.35) 



= -2(n - l)Mg/;+2(Mg) - (n' + 6n - 7)4+2(Mg), (2.36) 



P3i = 0, p32 = -4(Mg), (2.37) 



P33 = ^^^/n-2(Mg), P34 = |^^/n+2(Mg), (2.38) 

(n — zj (n + 2) 



P41 = ^In{Mq), P42 = -24(Mg), P43 = -4_2(Mg), P44 = -4+2(^5). (2.39) 
Transverse-traceless tensor modes /c„(t) yield, by virtue of Eqs. (2.10) and (2.19), 

Dn{Mq) = In{Mq), Vn > 3, (2.40) 

with degeneracy 2(n^ — 4). 

As far as ghost modes arc concerned, the decoupled mode /i(r) = ^/2(Mr) vanishes at 
the three-sphere boundary and hence yields 

D{Mq)=h{Mq), (2.41) 

with degeneracy 1, while scalar and vector ghost modes lead to 

Dr,{Mq)^I^^^{Mq), Vn > 2, (2.42) 

and 

D„(Mg) = 4(Mg), Vn > 2, (2.43) 

respectively, with degeneracy n^ for Eq. (2.42) and 2(n^ — 1) for Eq. (2.43). 
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Our D{Mq) equations can be re-expressed in a very helpful way by using repeatedly the 
identities for modified Bessel functions and their derivatives in the Appendix. Hence we 
find, on setting w = Mq, that 



D{w) 



(2.44) 



for the decoupled vector mode in Eq. (2.20), while coupled vector modes in Eq. (2.21) yield 



(n + 2) 
4n 



Moreover, the scalar modes ai,ei ruled by Eq. (2.26) yield 

^ = /.w(4W + l/.w), 

and the scalar modes 02, 62, 62 ruled by Eq. (2.30) lead to 

D(w) ^ . , ^ , , , , 5 ^ , , 
■-h{w)l3{w){l'^{w) + -h{w) 



6w^ 



w 



(2.45) 



(2.46) 



(2.47) 



III. MATRIX FOR COUPLED SCALAR MODES 

The hardest part of our analysis is the investigation of Eq. (2.31). For this purpose, we 
first exploit the formulae in Appendix B to find 

PU = In{w)-Wln{'w), P12 = Wl'^{w), P13 = wl'^{w) + nln{w) , pu = wl'^{w)-nln{w), (3.1) 



P21 = -{n^ - l)4(w), P22 = 2{wl'^{w) + 3/„(w)), 



(3.2) 



P23 = (n + 1) 



3(n + l) + 



2n(n- l)(n + 3) 



In{w) + 2 



w + 



{n- l)(n + 3) 



w 



(3.3) 



P24 = (n - 1) 



■ , 2n(n + l)(n-3) 
- 1) + —2 



In{w) - 2 



w + 



(n+ l)(n-3) 



w 



(3.4) 



P31 = 0, P32 = -In{w), 



(3.5) 
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P33 



(n+ 1) 
(n-2) 



1 + ^ . ^ 1 /n(«^) + 



w 



(3.6) 



P34 



jn-l) 
(n + 2) 



2n(n + l)^. 
1 H n I 



2(n + l] 



w 



I'nH 



(3.7) 



P41 = 34(w), P42 = -2In{w), 



(3.8) 



/ 2n(n- 1)\ ^ , , 2(n- 1) ^, 



w 



(3.9) 



P44 = - ( 1 + 4(^.) + '-^ILH. (3.10) 

The resulting determinant, despite its cumbersome expression, can be studied by introducing 
the variable 



y = 



which leads to 



where 



and hence 



Dn{w) = ^ ^ In{w){y - yi){y - y2){y - yi){y - yi), 



n n n w n w 

yi = , 2/2 = -, y3 = IT, 1/4 = — —, 

w w w 2 w 2 



^^^D.(.) ^ + (/»-^4(^)) 



(3.11) 



(3.12) 



(3.13) 



X {i» + (f + ^) inH) (^:h + (f - ^) ^'^(^)) • (3-14) 



IV. FIRST PAIR OF SCALAR-MODE (-FUNCTIONS 



Equations (2.40)-(2.47) and (3.14) are sufficient to obtain an integral representation of 
the C-function, the residues of which yield all heat-kernel coefficients. This topic is described 
in great detail in the existing literature (see, for example, Refs. [11] and [15]) and hence we 
limit ourselves to a very brief outline before presenting our results. 



12 



Given the elliptic operator P acting on physical fields defined on the m-dimensional 
Riemannian manifold A4, one can build the associated heat kernel U{x,y;t) and the cor- 
responding integrated heat kernel (bundle indices are not written down explicitly, but the 
fibre trace tr takes them into account) 

Tri2e-*-f'= / tTU{x,x;t)^ (Tx, (4.1) 
Jm 

which has the asymptotic expansion, as t — > 0"*", 

oo 

Tri2e-*^ ~ (47rt)-^ J2 ^^t^ ■ (4-2) 

fe=0 

The Ak coefficients are said to describe the global asymptotics in that they are obtained by 
integration over A4 and its boundary B of local geometric invariants built from the Riemann 
curvature of Ai, gauge curvature, extrinsic curvature of B, potential terms in P and in the 
boundary operator expressing the boundary conditions. On the other hand, since the 
function of P is related to the integrated heat kernel by an inverse Mellin transform [1,3,11]: 

Cp(s) = Trr^.p-^ ^ ^ = -i- f-' (Tri.e-*^)dt, (4.3) 

the global heat-kernel coefficients in the asymptotic expansion (4.2) can be also obtained 
from the residues of Cp{s) [15]. 

Moreover, since the function occurring in the equation obeyed by the eigenvalues by 
virtue of the boundary conditions admits a canonical-product representation [1,3], one can 
also express Cp{s) as a contour integral which is eventually rotated to the imaginary axis. 
The residues of the latter integral yield therefore the Ak coefficients used in evaluating 

2 

one-loop effective action and one-loop divergences. 

In our problem the P operator is the Laplacian on the Euclidean four-ball acting on 
metric perturbations. Equations (2.40)-(2.47) correspond to a famihar mixture of Dirichlet 
and Robin boundary conditions for which integral representation of the C-function and heat- 
kernel coefficients are immediately obtained. New features arise instead from Eq. (3.14), 
that gives rise, at first sight, to four different (^-functions. On studying the first line of Eq. 
(3.14), we can exploit the work in Ref. [16] and the uniform asymptotic expansion of Bessel 
functions and their first derivatives (see Appendix B) to say that the integral representation 
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of the resulting (^-function reads as 

(±(5) ^(^^y^-i2s-2) r z-'^^\og \z-P^^^'>(znI'^{zn)±nUzn)]] . (4.4) 
~, Jo dz I \ /J 



n=3 



With our notation, (3+{n) = n, j3^{n) = n + 2, where these factors are fixed by the leading 



behaviour of the eigenvalue condition as z — > 2; the uniform asymptotic expansion of 
modified Bessel functions and their first derivatives (see Appendix B) can be used to find 
(hereafter r = t{z) = (1 + from Eq. (B8)) 



00 



where 

Pfc,±(r) = (1 ± r)-i (t;fc(r) ± ^^(r)) , (4.6) 

for all A; > 1, and 

V k=l J k=l 

Thus, the ^-functions (4.4) obtain, from the first pair of round brackets in Eq. (4.5), the 
contributions (cf. Ref. jl^ ) 

A^s) ^ .'^-^log ' + ^'^'^ j , (4.9) 

where in the denominator of the argument of the log arises, in Eq. (4.9), from the extra 
z~'^ in the prefactor the definition (4.4). Moreover, the second pair of round 

brackets in Eq. (4.5) contributes Yl'jLiAj^±{s), having defined 

A,As) - f;n-(2^+^-^)i^ Hdz .-^^^T,,i(r(.)), (4.10) 

n=3 ^ 

where, from the formulae 

ri,±=Pi,±, (4.11) 



T2,±=P2,±-\pU, (4.12) 
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T3,± = P3,± - Pi,±P2,± + IpI±, (4.13) 

we find 

T,^^-lr±lr'-^T\ (4.14) 

' 8 2 24 ^ ^ 

T,^ = -Ir^ ± + ir^ ^ ^r' + -r^ (4.15) 
^'^ 16 8 8 8 16 ' ^ ^ 

21 3 , 3 4 509 . 25 6 21 7 , 15 8 1105 « 

Ta ± = ± -r^ + t —r^ + ± — r^, (4.16) 

128 8 640 12 128 8 1152 ' ^ ^ 

and hence, in general, 

ri,±(r)=i;/i^'^v«. (4.17) 

a=j 

We therefore find, from the first line of Eq. (3.14), contributions to the generahzed 
(■-function, from terms in round brackets in Eq. (4.5), equal to 

oo 

5C^(s) =a;o(5)Fo±(.) + J]a;,(.)if (s), (4.18) 

where, for all \ — 0,j, 

oo 

uxis) = J]n-('^+^-') = Ch{2s + a - 2; 3) = CrC^s + A - 2) - 1 - 2-^^'+^-^\ (4.19) 

n=3 

while, from Eqs. (4.8)-(4.10), 

F+{s)=^^^ rdzz-'^-^\og(l + {l + zY'A, (4.20) 

71 Jq OZ \ / 

F,-is) ^ -2i^ 2 - ^ois) = -1 - i^o^(«), (4.21) 

F±^,^^i^I^f2LHs.a,0)fi-^\ (4.22) 

having set (this general definition will prove useful later, and arises from a more general 
case, where r" is divided by the 6-th power of (1 ± r) in Eq. (4.17)) 

L^{s, &) = _ ^ys(^^ ^ ^ys b{l ± t)-'''^ - ar-\l ± r)-'')dr. (4.23) 
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Moreover, on considering 



and changing variable from z to r therein, all L-type integrals above can be obtained from 

Q{a,f3,-f)= [ T^^il-Tfil + rydr. (4.25) 
Jo 

In particular, we will need 

LUs) = -Q{2s,-s,-s-1), (4.26) 



L+(s, a, b) = bQ{2s + a, -s, -s-b-1) - aQ{2s + a - 1, -s, -s - 6), (4.27) 
where, from the integral representation of the hypergeometric function, one has |17|] 

g(a, /3, 7) = r(« + l)r(/3 + l) ^^_^^ a + l;a + P + 2;-l). (4.28) 



T{a + p + 2) 



For example, explicitly. 



Ltis) = - ^^^'rJ+2) + 1, 2. + 1; . + 2; -1). 

Now we exploit Eqs. (4.4), (4.5) and (4.18) to write 

^-(2s-l) 



(4.29) 



Clis) = SCXis) + 



SmTTS 



IT 



71=3 



2{l + z^) 



Hence we find 



C(0) = lim 



(4.30) 



(4.31) 



The first limit in Eq. (4.31) is immediately obtained by noting that 



limL+(.) = -log(2), 

s— >o 



and hence 



lima;o(s)FQ'~(s) = lim 



C„(2.-2;3)(!i!ii^tt(.) 



TT 



(4.32) 



(4.33) 



To evaluate the second limit in Eq. (4.31), we use 



lim L'''(s, a, 0) = — 1, 



(4.34) 
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and bear in mind that ujj{s) is a meromorphic function with first-order pole, as s — > 0, only 
aX j = 3 by virtue of the limit 



lim 



7- 



(4.35) 



Hence we find (see coefficients in Eq. (4.16)) 

OO / . N oo 



3i 



9 ^ ^ a 



(3,+) 



a=3 



720' 



(4.36) 



while, from Eqs. (4.30) and (4.28), 



lim(ci(.) - SC(s)) = lim (\U2S - 2; 3) + -^£g^C«(2. - 3; 3) 



.0 \ 4 

_ _5 1079 
~ ~4 ^ 240 ■ 

From Eqs. (4.31)-(4.37) we therefore obtain 

5 1079 1 



(4.37) 



(3,+) 



a=3 



146 
~A5' 



(4.38) 



An analogous procedure leads to 

/- 5 1079 



2^''« 90' 



(4.39) 



a=3 



where the +5 contribution results from Eq. (4.21) when exploited in Eq. (4.18). These 
results have been double-checked by using also the powerful analytic technique in Ref . 1^ ■ 



V. FURTHER SPECTRAL ASYMPTOTICS: ELLIPTIC AND NON-ELLIPTIC 
PARTS 

As a next step, the second line of Eq. (3.14) suggests considering (^-functions having the 
integral representation (cf. Eq. (4.4)) 

n=3 

/■OO Q 
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± n IJzn] 



(5.1) 



To begin, we exploit again the uniform asymptotic expansion of modified Bessel functions 
and their first derivatives to find (cf. Eq. (4.5)) 



where we have (bearing in mind that uq — vq — 1) 

2t 
(l~r2) 

for all A; > 1. Hereafter we set 



'^fe,±(T) = Ukir) + _^ (^{vk-iir) ± TUk-iir)^ , (5.3) 



k=l 



and rely upon the formula 

log(l + Q)~^(-l)^+Y (5.5) 



k 

k=l 



to evaluate the uniform asymptotic expansion (cf. Eq. (4.7)) 



w / ^ — ' w 

k=l / k=l 



The formulae yielding Rk,± from rk,± are exactly as in Eqs. (4.11)-(4.13), with T replaced 
by R and p replaced by r (see, however, comments below Eq. (5.10)). Hence we find, bearing 
in mind Eq. (5.3), 



' ^ ^ V 8 8 24 24 y 



(5.7) 



^ ^ ^ * 16 8 16 4 16 8 16 ' ' ^ ^ 

,3/1721 3 441 4 597 5 1033 g 239 , 

^ ^ V 384 128 320 960 80 

28 2431 221 _ 1105 1105 ,\ 

^ 5 ^ 576 192 384 1152 J' ^ ^ 

and therefore ^ 

i2,,±(r(^)) = (lTr)-^f;C^'±V», (5.10) 

a=j 

where, unlike what happens for the Tj^± polynomials, the exponent of (l^r) never vanishes. 
Note that, at r = 1 (i.e. ^ = 0), our rk,+{T) and Rk^+ir) are singular. Such a behaviour 
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is not seen for any of the strongly elliptic boundary-value problems (see third item in Ref. 



ll|). This technical difficulty motivates our efforts below and is interpreted by us as a clear 
indication of the lack of strong ellipticity proved, on general ground, in Ref. p. 

The Cb{s) function is more easily dealt with. It indeed receives contributions from terms 
in round brackets in Eq. (5.2) equal to (cf. Eq. (4.9) and bear in mind that /3_ — /9+ = 2 in 
Eq. (5.1)) 



B_is) ^ y,-(2.-2)(^ r^z.-^logf^ 
tri Jo dz \ 



. Z2 

n=3 

fsinvrs) /""^ . _2, <9 , 1 1 



Ms)- dz z log = --ujois), (5.11) 







dz 



and Yl'jLi Bj-{s), having defined (cf. Eq. (4.10)) 

BjAs) = ^,(g) ^ / dzz 2^— i?^._(r(z)) = ujis) ' ^ ^ 2^L+(g,a,j)Cp- ). 

a=j 

(5.12) 

On using the same method as in Sec. IV, the formulae (5.1)-(5.12) lead to (we find 
L+(0,a,3) = — |, independent of a, below) 

a=3 

a result which agrees with a derivation of Cij(O) relying upon the method of Ref. |l4| . 

Although we have stressed after Eq. (5.10) the problems with the Csi^) part, for the 
moment let us proceed formally in the same way as above. Thus we define, in analogy to 
Eq. (5.11), 

BAs)^Ms)^^^ r dzz-'^^ log (^-riz)) , (5.14) 

TT Jo dz \t{z) J 

and, in analogy to Eq. (5.12), 

i?,,+ (.) ^ io.is)^^^ r dz z-'^^R,Ar{z)). (5.15) 
TT Jq dz 

In order to make the presentation as transparent as possible, we write out the derivatives of 
Rj,+- On changing integration variable from z to r we define 

d_ 

and we find the following results: 



Q(r) = Tr^i,+ M, (5-16) 



C.{r) = (1 - r)- C^-\r- + \r' -\r' ] . (5.17) 
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G.'^(t) = (1-r) \ -r^T — ^ -r^T — —T + 



/1721 2 
( 128 ^ 


441 o 
- 32" 


+ 


1635 4 
128 ^ 


163 , 


221 „ 


5083 10 


+ 


1105 11 


1105 


16" 


128 


r - 

32 


128 



+ 4223^8 _ _ ^^10 ^ ^^11 _ ^^12 J ^ ^^^g^ 



so that the general expression of Cjij) reads as 

4i ^(i) a 

^^•(^) = E n-rV+i ' Vj = 1, . . oo . (5.20) 



(1 - 

These formulae engender a Cs (0) which can be defined, after change of variable from z to r, 
by splitting the integral with respect to r, in the integral representation of Csl'^)' according 
to the identity 

/•I /"M /"i 

/ dr = dT+ dr, 
Jo Jo J n 

and taking the limit as ^ 1 a/i(er having evaluated the integral. More precisely, since the 
integral on the left-hand side is independent of |U, we can choose /i small on the right-hand 
side so that, in the interval [0,^] (and only there!), we can use the uniform asymptotic 
expansion of the integrand where the negative powers of (1 — r) are harmless. Moreover, 
independence of also implies that, after having evaluated the integrals on the right-hand 
side, we can take the fj, ^ 1 limit. Within this framework, the limit as — > 1 of the second 
integral on the right-hand side yields vanishing contribution to the asymptotic expansion of 

With this caveat, on defining (cf. Eq. (4.25)) 



g^(a,/9,7)=/ T^^il-rni + rydT, (5.21) 







we obtain the representations 
(sin ns) 



( sm Tis I r "1 
B+{s) = -uoisy - Q^{2s, -s - 1, -s) + Q^(2s, -s, -s - 1) - Q^(2s - 1, -s, -s) , 

TT L J 

(5.22) 



sm TTS 



B,,^{s) = -uj.is)^- L 5^ Ki^^Q,{2s + a,-s-j-l,-s). (5.23) 



TT 

a=j-l 
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The relevant properties of Q^(a;,/3,7) can be obtained by observing that this function is 
nothing but a hypergeometric function of two variables [l7|, i.e. 

/3, 7) = -^Fi(a + 1, -(3, -7, a + 2; Ai, -fi). (5.24) 
a + 1 

In detail, a summary of results needed to consider the limiting behaviour of Csi^) as s ^ 



IS 



Ms)- -5s + 0{s^), (5.25 

TT 



, , (sin ns) 1 ^ - ^ , / n 
^ ~ 3 + + 0(^2), (5.26) 

TV 2 



lim Q^(2s, -s - 1, -s) ~ -- + 0(s°), (5.27) 



lim Q^(2s, -s, -s - 1) ~ log(2) + 0(s), (5.28) 

limQ^(2s-l,-s,-s) ~ -^ + 0(s), (5.29) 
M^i 2s 

limQ^(2s + a,-s-j-l,-s) = ^^^^^=-^E^^^^^^F(a + 2s + 1, s, a - j + s + 1; -1) 

~ ^^^ + 6,,o(«) + 0(s), (5.30) 

S 

where 

6,, 1 = -1 - 22-^' + CnU - 2)(1 - ^,,3) + lS,,s, (5.31) 



fe,,„i(a) = ^ . /,, (! - '^aj-i), (5.32) 



( -1)^'+^ r(a + l) 
j! r(a-j + l) 

and we only strictly need 63,0 («) which, unlike the elliptic cases studied earlier, now depends 
explicitly on a. In our case, we find 

12 67 95 

&3,o(2) = -, 63,0(3) = --log(2), 63^0(4) = --41og(2), 63,0(5) = - -10 log(2), 

14-3 139 2433 

^>3,o(6) = -f- 201og(2), &3,o(7) = — - 351og(2), 63,0(8) = ^ " 561og(2), 

1971 8429 
^^3,0(9) = ^ - 841og(2), 63,0(10) = ^ - 1201og(2), 
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&3.o(ll) = ^ - 1651og(2), 63,0(12) = - 2201og(2). 

56 504 

Remarkably, the coefficient of ^ in the small-s behaviour of the generahzed (^-function 
Cb(s) is zero because it is equal to 

12 ^ 12 

5^63,_i(a)irf = -5^a(a- l)(a-2)irf , 

a=2 a=3 

which vanishes by virtue of the rather peculiar general property 

E rfa-t+l) ^"' = - = 0, V, = 1, . . ., 00, (5.33) 

a=j ^ •' ' a=j 1=0 

and hence we find eventually 

"-v I I I / V J "-X I _ 

j=l a=j-l 



5 

^4 ' 


1079 
^ 240 


5 

+ 2 


1 

a=2 


5 

4 + 


1079 
240 ^ 


599 
720 


296 
45 ' 



(5.34) 

because the infinite sum on the ffist line of Eq. (5.34) vanishes by virtue of Eqs. (5.32) and 
(5.33), and exact cancellation of log(2) terms is found to occur. 

To cross-check our analysis we remark that, on applying the technique of Ref. one 
finds 

CJ(0, = -15 + H2_^ = 5!, (5.35) 

^ 4 240 720 90 ^ ^ 

where results from working in the 00 and w ^ limit in 



{I'nM + ( 



/'w n 
2^w 



IJw] 



on the second line of Eq. (3.14); such a term then reduces to [l'^{w) + ^/„(t(;)). A possible 
interpretation of the discrepancy between (5.34) and (5.35) is that, when strong ellipticity 
is violated, prescriptions for defining a C(0) value exist but are inequivalent. 

Remaining contributions to C(0)) being obtained from strongly elliptic sectors of the 
boundary- value problem, are instead found to agree with the results in Ref. 0, i.e. 

278 

C(0) [transverse traceless modes] = , (5.36) 

45 



494 

C(0) [coupled vector modes] = (5.37) 
22 



C(0) [decoupled vector mode] = —-pr, 



(5.38) 



C(0)[scalar modes(ai, Ci; 02, 62, 62)] = —17, 



(5.39) 



C(0) [scalar ghost modes] 



149 



(5.40) 
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C(0) [vector ghost modes] 



77 

90' 



(5.41) 



C(0) [decoupled ghost mode] 



5 



(5.42) 
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VI. CONCLUDING REMARKS 

We have obtained the analytically continued eigenvalue conditions for metric perturba- 
tions on the Euclidean four-ball, in the presence of boundary conditions completely invariant 
under infinitesimal diffeomorphisms in the de Bonder gauge and with a parameter set to 1 
in Eq. (1.10). Second, this has made it possible to prove, for the first time in the literature, 
that only one sector of the scalar-mode determinant is responsible for lack of strong ellip- 
ticity of the boundary- value problem (see second line of Eq. (3.14) and the analysis in Sees. 
IV and V). The first novelty with respect to the work in Ref. |6| is a better understanding of 
the elliptic and non-elliptic sectors of spectral asymptotics for Euclidean quantum gravity. 
Moreover, as far as we know, the detailed spectral asymptotics for ^-functions of Sees. IV 
and V was missing in the literature. We have also shown that one can indeed obtain a 
regular C-function asymptotics at small s in the non-elliptic case by virtue of the remarkable 
identity (5.33), here obtained for the first time. Our prescription for the C(0) value differs 
:'rom the result first obtained in Ref. 0, where, however, neither the strong ellipticity issue 
^ nor the non-standard spectral asymptotics of our Sec. V had been considered. 

From the point of view of general formalism of Euclidean quantum gravity, three alter- 
native pictures seem therefore to emerge: 

(i) The remarkable factorization of eigenvalue conditions, with resulting isolation of elliptic 
part of spectral asymptotics (transverse-traceless, vector and ghost modes, all modes in 
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finite-dimensional sub-spaces and three of tlie four equations for scalar modes), suggests 
trying to re-assess functional integrals on manifolds with boundary, with the hope of being 
able to obtain unique results from the non-elliptic contribution. If this cannot be achieved, 
the two alternatives below should be considered again. 

(ii) Luckock boundary conditions which engender BRST-invariant amplitudes but are 
not diffeomorphism invariant P]. They have already been applied by Moss and Poletti jigj ]. 

(iii) Non-local boundary conditions that lead to surface states in quantum cosmology and 



pseudo-differential operators on metric and ghost modes [2l|. Surface states are particularly 
interesting since they describe a transition from quantum to classical regime in cosmology 
entirely ruled by the strong ellipticity requirement, while pseudo-differential operators are a 
source of technical complications. 

There is therefore encouraging evidence in favour of Euclidean quantum gravity being 
able to drive further developments in quantum field theory, quantum cosmology and spectral 
asymptotics (see early mathematical papers in Refs. [2^, [2^) in the years to come. 

APPENDIX A: STRONG ELLIPTICITY 

For an operator of Laplace type, the boundary-value problem is strongly elliptic with 
respect to the cone C — R+ if, for any cotangent vector u on the boundary B, for any 
A G C — R+, for any pair (u, A) 7^ (0, 0), there exists a unique solution tp of the differential 
equation (r being the geodesic distance to the boundary B) 

Q2 



Qj.2 

subject to the asymptotic condition 



+ UkU - A 



ip{r) = 0, (Al) 



lim ¥?(r) = 0, (A2) 



and to the boundary conditions (here ^{r) = '^^ with a = yuj^vF—X) 

7i^{r = 0) = i/jo, iT(f{r = 0) + (/ - n)<p'{r = 0) = i^i, (A3) 

where vr is the same projector as in Eq. (1.1), iT is the leading symbol of that part of 
the boundary operator which involves tangential derivatives, while i/jq and '01 are arbitrary 
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boundary data. Eventually, all this is equivalent to proving positivity of the matrix I\/ukU^ — 
iT 0. 



APPENDIX B: BESSEL FUNCTIONS 

In Sec. II we exploit the following identities obeyed by modified Bessel functions of first 
kind: 

Tl 

4+iH = /;H--/„H, (bi) 

w 



Tl 

In-liw) = I'M + -InM, (B2) 
W 

, , / 2n(n + l)\ , , 2(n+l) , , , 

\ J w 

r / N f 2n(n-l)\^, , 2(n-l)^, , , 
\ J w 

, , , 2(n + l) / n(n + 2)\ , , / 2(n+l)(n + 2)\ , , , 

, , , 2(n-l) / n(n-2)\ , , ( 2(n-l)(n-2)\ , , , 

In Sees. IV and V we use the uniform asymptotic expansion of modified Bessel functions 
4 first found by Olver 



iJzv) 



^2^(1 + 

where 

T={1 + z'yl, ri = {l + z'Y^ + log ( ^j== ) . (B8) 



,1 + v^rT^y ■ 

This holds for 1/ — >■ 00 at fixed z. The polynomials Mfc(r) can be found from the recurrence 
relation 

nk+i{r) = ]-t\1 - t^)u'^{t) + 1 [ dp (1 - bp'^Wip), (B9) 
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starting with Uo{t) = 1. Moreover, the first derivative of I^, has the following uniform 
asymptotic expansion at large u and fixed z: 



e^^ {l + z^ 



VkVr) 



k=l 



with the Vk polynomials determined from the Uk according to [Ji 

'1 



Vk{T) = Mfc(r) + r(r^ - 1) 
starting with Vo{t) = Uq{t) = 1. 



(BIO) 



(Bll) 
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